Single- and double-vortex vector solitons in self-focusing nonlinear media 
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We study two-component spatial optical solitons carrying an angular momentum and propagating 
in a self-focusing saturable nonlinear medium. When one of the components is small, such vector 
solitons can be viewed as a self-trapped vortex beam that guides either the fundamental or first-order 
guided mode, and they are classified as single- and double-vortex vector solitons. For such composite 
vortex beams, we demonstrate that a large-amplitude guided mode can stabilize the ring-like vortex 
beam which usually decays due to azimuthal modulational instability. We identify different types of 
these vector vortex solitons and demonstrate both vortex bistability and mutual stabilization effect. 
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I. INTRODUCTION 

Vortices are fundamental localized objects which ap- 
pear in many branches of physics Q . In fluid mechan- 
ics, coherent structures in the form of vorticity filaments 
are central dynamical objects to understand most fluid 
flows and particularly fluid turbulence. More recently, 
the study of vortices in Bose-Einstein condensates re- 
vealed many intriguing properties of superfluids created 
by ultra-cold atoms 2\. Different types of vortices can 
also be found and identified in optics; one of the simplest 
objects of this kind is a phase singularity in an optical 
wave front which is associated with aphase dislocation 
carried by a diffracting optical beam 0, Q • 

In self-focusing saturable nonlinear media, optical vor- 
tices can exist as self-trapped ring-like optical beams 
with zero intensity at the center carrying a phase sin- 
gularity @|. However, due to the self-focusing nature of 
nonlinearity such ring-like vortex beams become unstable 
to azimuthal perturbations, and they decay into several 
fundamental optical solitons flying off the main ring Q. 
This effect was observed experimentally in different non- 
linear media, including the saturable Kerr-like nonlin- 
ear media, biased photorefractive crystals, and quadratic 
nonlinear crystals operating in the self-focusing regime 
(see details and references in Ref. Q). There are known 
several ways to stabilize this azimuthal modulational in- 
stability, including the vortex stabilization in the pres- 
ence of a large-amplitude beam guided by it and the 
stabilizating effect of partial incoherence of light on the 
vortex |3 ■ 

When a self-trapped vortex beam guides a large- 
amplitude fundamental beam, it creates together with 
the guided beam a composite object in the form of a vec- 
tor soliton 0- Mutual coupling between the fundamental 
beam and the vortex-carrying beam can create different 
novel types of composite vector solitons carrying an angu- 
lar momentum [lot fill fl^ | ■ The properties of such vector 
vortex solitons can differ substantially from the proper- 
ties of one-component scalar vortices and scalar solitons. 
In particular, the mutual coupling between the beams 
can modify dramatically the vortex properties and, in 



particular, can suppress the development of azimuthal 
instability [8|. 

In this paper, we analyze the existence, general prop- 
erties, and stability of the vector vortex solitons in a 
self-focusing saturable nonlinear medium. We study two 
types of such two-component composite vortex solitons, 
i.e. single-vortex vector solitons that can be considered as 
the vortex-induced waveguide that guides a fundamental 
mode, and double-vortex vector solitons, when the local- 
ized field is similar to the first-order guided mode, be- 
ing a vortex beam by itself. For some of the cases, we 
demonstrate that a mutual incoherent coupling between 
the vortex waveguide and a large-amplitude guided mode 
it guides can provide a strong stabilizing mechanism for 
stable or quasi-stable two-component vortex solitons to 
exist in such media, in agreement with the recent obser- 
vation 13] of the stabilizing mechanism of the mutual 
coupling between different components of the composite 
vortex beam. 



II. MODEL 

In order to study the vector vortex solitons, we con- 
sider the interaction of two mutually incoherent optical 
beams propagating in a self-focusing nonlinear saturable 
medium. The evolution equations for two incoherently 
interacting beams can be presented in the following di- 
mensionless form, 



du 

i— + A ± u - 

az 

dv 

i— + Aj_v - 

az 



(\u\ 2 + n\v\ 2 )u 
1 + a{\u\ 2 + M 2 ) 

(H 2 + n\u\ 2 )v 
I + a(M 2 + M 2 ) 



= o, 

= 0, 



(1) 



where u and v are the dimensionless amplitudes of the 
fields, the parameter a characterizes the nonlinearity sat- 
uration effect, and the mutually incoherent interaction 
between the modes is described by the coupling param- 
eter ji. The spatial coordinate z is the beam propaga- 
tion direction, and Aj_ stands for the transversal part 
of the Laplace operator in the cylindrical coordinates, 



2 



defined through the radial, r = (x 2 + y 2 ) 1 / 2 , and angu- 
lar, <p = tan _1 (y/x), coordinates. Model Q provides a 
straightforward generalization to a number of important 
cases studied earlier. In particular, the limit a — > cor- 
responds to the Kerr medium with cubic nonlinearity [T^ 
where all self-trapped beams may undergo collapse insta- 
bility. The case of the saturable nonlinearity at fi = 1 
corresponds to the incoherent beam interaction in pho- 
torefractive nonlinear media [lot fill fl4| . We look for 
stationary solutions of the system that describe a ra- 
dially symmetric single-charged vortex beam in the field 



u(r, <p;z) 



u(r)e tlf >e iz , 



(2) 



where the amplitude u(r) vanishes for r — > oo. We as- 
sume that the vortex (JSJ) guides (or is coupled to) the 
second beam, 



v(r, (j); z) = v(r)e 



fc- n 



(3) 



where f3 is a dimensionless ratio of the propagation con- 
stants. In Eqs. (0), (J3J) the functions u{r) and v(r) are the 
radial envelopes of the interacting fields, and I (I = 0, ±1) 
is the angular momentum of the guided mode. Equations 
for the stationary envelopes are given by 



A r U 



I 2 

A r v -v 



(u 2 + fiv 2 )u 
1 + cr(u 2 + v 2 ) 
(v 2 + (M 2 )v 



0. 

= 0, 



(4) 



r* 1 + o~(u 2 + v 2 ) 
where A r is the radial part of the Laplace operator, 



A - 1 d ( d 
r r dr \ dr 



The radially symmetric, spatially localized solutions of 
the system Q describe different types of two-component 
composite solitons carrying an angular momentum, and 
they can form either single- or double- vortex vector soli- 
tons. In a two-dimensional geometry, such solutions can 
only be found numerically. 



III. SINGLE- VORTEX VECTOR SOLITONS 

We consider the self-trapped vortex beam created in 
the main field u, with the asymptotic behavior u{r) — ► 
as r — * oo. Besides that, we require that the incoherently 
interacting component v(r) describes a localized mode, 
i.e. v(r) — » and therefore (3 > 0. At the origin (r = 0), 
the boundary condition for the vortex is u — and, if 
we seek the single-vortex solution with the fundamental 
guided mode (I — 0), the corresponding boundary condi- 
tion for v(r) is dv/dr — (or v = vq). 

We find localized solutions numerically, by means of 
the relaxation technique. In Fig. ^ we show some ex- 
amples of two-component localized solutions which de- 
scribe a fundamental (no nodes) beam, guided by the self- 
trapped vortex that create together a single-vortex vector 
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FIG. 1: Examples of single- vortex vector solitons created by 
the vortex beam it and the fundamental mode v it guides. 
Labels correspond to the points marked in Fig. |5] For points 
B and D, we show two different solutions named with roman 
numbers I and II, that exist due to bistability. The model 
parameters are: A (/i = 1.10,/? = 1.38), B (/i = 1.10, /3 = 
1.44), C (fi = 1.17, p = 1.38), and D (// = 1.17, (3 = 1.47). 



soliton. The existence domain for such solutions has been 
calculated numerically for a special case a = 0.5, and it is 
shown in Fig. on the parameter plane (/?, /x). The exis- 
tence region for single- vortex vector solitons is composed 
of two regions which overlap in a triangular-shaped do- 
main shaded with a different intensity in Fig. ^ In this 
intersection, two types of vortex solutions mark the fa- 
miliar bistability phenomenon. The existence domain is 
restricted by the solid curves which describe some spe- 
cific cutoff boundaries. Close to the lower cutoff, one of 
the components becomes small: the fundamental field in 
the left region (see the cases A and D-I), and the vortex 
in the right region (cases B-II and C). The other compo- 
nent is therefore only weakly distorted. However, when 
the propagation constant is close to the upper cutoff 
(cases B-I and D-II) , the amplitudes of both components 
become comparable, affecting strongly each other. For 
small values of the coupling parameter /i, the existence 
region is limited from above by the constant value 1/cr 
(in our example, 1/cr = 2), the value which can be easily 
explained by a simple qualitative analysis. Indeed, ac- 
cording to Eqs. 10}, the existence of bounded stationary 
states for both u and v components requires the following 
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FIG. 2: Existence domain of the single-vortex vector soli- 
tons on the parameter plane (/3, /i) at a — 0.5. Inset shows 
an enlarged region of bistability where two different types of 
vortex-mode solutions coexist. 
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FIG. 3: Bifurcation diagram for the two-component vector 
vortex solitons, shown together with the families of the scalar 
fundamental and vortex solitons at /i = 1.1. Thin curves 
represent the powers P u and P v of each of the scalar solitons; 
the thick curve is the total power Ptotai of a composite vortex 
soliton. Points Oi and O2 are the bifurcation points where 
the vector soliton emerges, and the intersection points B-I 
and B-II correspond to the example of bistable single-vortex 
solutions for the point B {fi = 1.44) shown in Fig. 

conditions to be valid, 

f u 2 + fiv 2 } 
max < — — > > 1, (5) 

1 + a{u 2 + v 2 ) J w 

f v 2 + ^u 2 \ 

max < — — > > p. (b) 

[1 + a(u 2 + v 2 ) J 

Considering /i — > 0, we obtain from Eq. @) that the 
existence of localized solution requires that 

r v 2 ) 

P < max <^ 7T- > . (7) 

[ 1 + a{u 2 + v 2 )j w 

The right-hand-side term vanishes when the amplitude of 
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FIG. 4: Examples of the vortex propagation dynamics in the 
bistability domain. Shown are the field intensity profiles as 
gray-scale images at several propagation distances. Top: the 
components of the B-I vector soliton. Bottom: the compo- 
nents of the B-II vector soliton (see Fig. 1 and Fig. 3. ). 



the component v vanishes at the cutoff, but it approaches 
the value 1/er when the amplitude of the field v becomes 
large. This analysis is valid for any type of the guided 
mode. 

In order to describe the bistable vector solitons, in 
Fig. we display the bifurcation diagram of the two- 
component vortex-mode localized solutions, for the par- 
tial and total beam powers. Composite vortex-mode soli- 
tons presented by the power dependence Ptotai originate 
at the bifurcation point Oi where the mode v is small 
and can be described by the linear theory. For larger 
value of (3, this curve bends, and then it merges with the 
other partial power curve P v at the bifurcation point O2 
(see Fig. 121 . The bistable solutions B-I and B-II, which 
are presented in Fig. 1 being related to point B in Fig. 2, 
correspond to a single value of the propagation constant 
/3 in the bistability domain. Importantly, two solutions 
have different stability properties, and only one of them 
is stable, as shown in Fig. In general, the solutions 
belonging to the left region in the domain are unstable 
and those belonging to the right zone are stable. 

The incoherent interaction between the vortex beam 
and the localized mode it guides has the character of at- 
traction, and it may provide an effective physical mech- 
anism for stabilizing the vortex beam in a self-focusing 
nonlinear medium. Indeed, a scalar vortex beam becomes 
unstable in a self-focusing nonlinear medium due to the 
azimuthal modulational instability. In this case, the vor- 
tex splits into the fundamental beams that fly off the 
main vortex ring. On the other hand, the bright soli- 
tons are known to be stable in such media. We expect 
that a mutual attraction of the components in a two- 
component beam may lead to a counter-balance of the 
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FIG. 5: Vortex stabilization due to mutual interaction. 
Shown are the field intensity profiles as gray-scale images at 
several propagation distances. Top: breakup of a scalar vor- 
tex soliton (v — 0). Bottom: both the vortex beam and 
the fundamental guided mode (approximated by a Gaussian 
beam) propagate together. Parameters are: a = 0.5, /i = 1, 
and (3 = 1.45. 



vortex instability by the bright component when its am- 
plitude is large enough. To confirm this idea, we consider 
a two-component composite structure consisting of a vor- 
tex beam together with the fundamental mode it guides, 
both described by Eqs. (JIJ at rj = 1 and fx = 1.0. To 
study the mode stability, we propagate the stationary 
soliton solutions numerically. In Fig. [5] we compare the 
vortex breakup for the scalar and vector systems. In the 
top row, we show the propagation of a vortex alone in 
the scalar model; the vortex breaks up into two solitons 
which fly away after some distance. In the bottom row, 
we show the propagation of two coupled components (the 
vortex and bright mode it guides) . Due to a strong inco- 
herent coupling between the modes, the propagation of 
the vortex is stabilized for some propagation distances, 
so that the vortex breakup is delayed dramatically, as 
shown in Fig. flower row), or even become completely 
stable, similar to the case shown in Fig. ^Jlow row). We 
confirm this stabilization mechanism by performing our 
study for a Gaussian input beam of the bright compo- 
nent instead of the exact stationary state, as would be 
easier realizing in experiment. 



IV. DOUBLE- VORTEX VECTOR SOLITONS 

To study novel types of vector vortex solitons, we con- 
sider the vortex beam in the field u coupled to the first- 
order guided mode described by the solution with 
1 = 1, Thus, for the field v we look for a vortex- like lo- 
calized solution with the boundary condition v(r) — > at 
r = 0. Our analysis shows that there exist four different 
kinds of such solutions, and they form the families of the 
so-called double-vortex vector solitons. We show some 
examples of these solutions in Fig. and Fig. EI whereas 
the existence domains for all solutions of this type are 




FIG. 6: Examples of four types of double-vortex vector soli- 
tons calculated for a — 0.5, /j, = 0.3 and j3 = 1.5. All solutions 
correspond to the point A in Fig. Q 
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FIG. 7: Existence domain for different types of double- vortex 
solitons. Domain 1: only the solutions of type I; domain 2: 
only solution of type II; domain 3: solutions of types I, II and 
III; domain 4: solutions of types I, II and IV; domain 5: all 
four types of solutions exist. Two vertical thin dotted lines 
indicate the values of \x for which the bifurcation diagrams 
are presented in Figs. I§1 and 1101 



shown on the parameter plane (/?, \i) in Fig. 

In the region labelled with number 5, there exist four 
types of localized double- vortex solutions, as shown for 
the point A in Fig. [5] The first type (type I) is described 
by two rings where the ring in the field u is larger than 
the ring in the filed v; the opposite situation occurs for 
the solutions of type II (see Fig. EJ) . Other two types of 
solutions have one of the fields of a two-humped shape: 
either the field u (type III) or the field v (type IV). As 
the coupling parameter /i grows, the valley in the two- 
humped solutions becomes shallower, and it disappears 
when the solution crosses a thin dashed line on the exis- 
tence plane (Fig. 0), moving either to the domains 3 or 
4. In the upper domain (p > 1, domain 3), no solution of 
type IV can be found further up from this line, being con- 
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FIG. 8: Tristable solution found for \i — 0.5 and j3 = 1.5 
(point B in Fig.0 and the unique solution found for fj, = 0.7 
and (3 = 1.5 (point C). In all cases a — 0.5. 



sequently the boundary of its existence domain. Besides, 
the solution of type III degenerates in a single-humped 
solution and still exists up to the thin continuous line 
(see Fig. IU case B-III). The opposite occurs in the lower 
domain ((3 < 1, domain 4) , where the solitons of type III 
cannot be found to the right of the dashed line, and the 
solitons of type IV degenerate into single- humped modes. 

Single-humped solutions of types I and II both exist 
up to the continuous thin line, though the maxima of 
both fields u(r) and v(r) are shifted, approaching each 
other as the parameter \x increases (see the examples B-I 
and B-II in Fig. [SJl. Below this line, in the upper do- 
main (fi > 1, domain 3) the solutions II no longer exist, 
while the solutions I have a symmetric shape where the 
maxima of the fields approximately coincide (Fig. [HI case 
C). Dashed thick curve corresponds to the cutoff for such 
solutions, towards which the amplitude of the mode v be- 
comes small and does not influence the vortex mode in 
the field u. On the other hand, the continuous line is the 
upper cutoff where a vector soliton originates. 

Shape of these solutions can be explained qualitatively 
by a simple analysis. In fact, from Eq. (JSJ), assuming /i 
small and the maximum amplitude of v(r) much larger 
than that of u(r) (i.e. the parameter (3 close to the upper 
cutoff), we obtain the condition max[w 2 /(l + uv 2 )] > 1. 
Since v has a larger maximum, the former condition can- 
not be be satisfied unless the value of v is small at the 
position of the maximum of u, so that there still exists an 
effective potential with a relative maximum higher than 
1. That requires that both fields have maxima shifted 
enough to each other, explaining the shape of the solu- 
tions in region 5 of Fig. \7\ When fi is larger, however, 
this condition can be satisfied when the maxima of both 
solutions almost coincide. Because the field v becomes 
self-guided and the shifting is not possible (due to the 
coupling between both components), each field has to be 
localized in the region close to the other. This explains 
the behavior of the solutions when /i grows and the solu- 



FIG. 9: Bifurcation diagram for vector and scalar vortex 
solitons (at fj, — 0.7). P u and P v (thin lines) are the powers 
of scalar vortices created in each component separately, while 
-Ptotai (thick line) is the power of the vector soliton originated 
at the bifurcation points Oi and O2. 




propagation constant, (3 

FIG. 10: Bifurcation diagram for vector and scalar vortex 
solitons (at /1 = 0.53). P u and P v (dashed thin lines) are 
the powers of the scalar vortices created in each component 
separately. Other curves are labelled with a subscript corre- 
sponding to different types of vector solitons: type I (contin- 
uous thin line), type II (continuous thick line), and type III 
(dashed thick line). Points Oi and O2 are bifurcation points, 
and O3, O4 and O5 are critical points. An area near the 
bifurcation point O2 is shown in the inset. 



tion crosses the boundaries between region 5 and region 
3, and then moves to region 1. On the other hand, if 
is decreased, the component v{r) becomes smaller being 
guided by the component u(r). This explains the behav- 
ior of the solutions crossing the boundaries from region 
5 into region 4 and to region 2. 

Existence of different kinds of double-vortex vector 
solitons leads to multi-stability phenomena as well as 
more complicated bifurcation diagrams. In Fig. [5] and 
Fig. I1UI we present two examples of the bifurcation dia- 
grams for two different values of the coupling parameter. 
In the first case, at n = 0.7, only one solution for each 
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FIG. 11: Evolution of two types of stable double-vortex vec- 
tor solitons corresponding to point D (type II, top rows) and 
point E (type I, bottom rows). For each soliton, both com- 
ponents are shown for different values of the propagation dis- 
tance. For comparison, we show the corresponding unstable 
evolution when the component v is removed at the input. 



value of the propagation constant exists. At j3 = 1, both 
solutions of type I (j3 > 1, region 1 in Fig. EJ and type II 
{(3 < 1, region 2) merge together. In fact, for (3 = 1, it is 
deduced from Eqs. 10} that both solutions u(r) and v(r) 
become identical. 

For smaller values of /i, however, the existence of dif- 
ferent kinds of vortex solutions generate a variety of 
branches in the bifurcation diagram, as shown in Fig. 1101 
In this case, there exist both the bifurcation points Oi 
and O2 and the critical points O3, O4 and O5, and three 
of the four kinds of solutions exist for some values of (3. 
For other values of the coupling parameter different 
types of the bifurcation diagrams are obtained, and they 
all show a change of the solutions when a boundary be- 
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FIG. 12: Decay of the vortex vector solitons into 2+4 (top 
row) and 2+2 (bottom row) fundamental solitons. The upper 
case corresponds to the solution of type IV (point F; /i = 0.2, 
/3 = 1.41), and the lower case -to the solution of type I (point 
C; /i = 0.7, f3 = 1.5). 
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FIG. 13: Examples of the vortex instability scenarios. Top: 
decay of type IV vortex solitons corresponding to point G 
(fi = 0.49, /3 = 1.95). Bottom: quasi-stable propagation 
of the vortex with breathing components for type II vector 
solitons (point H; /j, = 0.5, (3 — 1.95). 
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tween different existence domains is crossed. 

In Fig. ^2 we show two examples of stable propaga- 
tion. The top rows present the evolution of the vector 
soliton corresponding to point D in Fig. \7\ For compari- 
son, we show unstable propagation of the first component 
alone, when the second vortex component is removed 
from the input. In the latter case, the vortex decays 
after propagating for some distance while the two com- 
ponent soliton remains virtually stable for much longer 
propagation distance owing to the coupling between both 
the components. In general, it is possible to achieve a 
double-vortex vector soliton which remains stable for an 
arbitrary distance provided we chose a state close enough 
to the upper cutoff, where one of the components has 
larger amplitude and the interaction between both com- 
ponents is strong. In the bottom rows, we show another 
example of the stable propagation corresponding to point 
E in Fig. □ 

Different types of the double-vortex vector solitons 
demonstrate a rich variety of the instability-induced sce- 
narios of their evolution. In Fig. 1121 we present two ex- 
amples of the vortex evolution with two characteristic 
scenarios of the vortex decay, producing either 2+4 (e.g. 
point F) or 2+2 (e.g. point C) fundamental solitons. Fi- 
nally, in Fig. El w e show two more complicated scenarios 
of the vortex instability to illustrate a variety of the pat- 
terns that can be observed. In the top row, where the 
input state corresponds to point G, the soliton decays in 
a complex way displaying a sequence of symmetric pat- 
terns. In the bottom row, where the input state corre- 
sponds to point H, the double-vortex soliton propagates 
in a quasi-stable way performing breathing radial oscil- 
lations. In all the cases discussed here, the topological 
charges of both the components are chosen with the same 
sign (+1). The case of the opposite charges has also been 
studied, and the similar evolution scenarios have been ob- 
served. 



V. CONCLUSIONS 

We have analyzed the existence and basic properties of 
the two-component composite optical beams carrying an 
angular momentum, the so-called vector vortex solitons. 
We have considered two major types of such solitons that 
propagate in self-focusing saturable nonlinear media and 
can be classified in the low-intensity limit through the 
fundamental and first-order localized modes guided by 
the main vortex beam. We have calculated the existence 
domains of such vortex composite solitons and studied 
numerically their stability to weak perturbations. In par- 
ticular, we have revealed a novel mechanism for stabiliz- 
ing the vortex azimuthal modulational instability by a 
co-propagating guided mode of a large amplitude. We 
have demonstrated, for the first time to our knowledge, 
the existence of bistable composite double- vortex solitons 
as well as studied their instability-induced dynamics. We 
believe that similar results can be obtained for compos- 
ite vortices in other types of nonlinear models describing 
the mutual coupling between several fields, and our re- 
sults can be useful for other fields such as the physics of 
the multi-species Bose- Einstein condensates of ultra-cold 
atoms. 
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